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Abstract. 

We present a comparison of the gas and dust properties of the dense interstellar matter in six nearby star— forming 
regions (d <500 pc): p Oph, Cha I, R CrA, IC 348, NGC 1333, and Orion. We measure from Chandra and 
XMM— Newton observations the X-ray absorption toward pre-main sequence stars (PMS) without accretion disks 
' (i.e., Class III sources) to obtain the total hydrogen column density Nh,k- For these sources we take from the 

, literature the corresponding dust extinction in the near— infrared, Aj, or when unavailable we derive it from SED 

fitting using the available DENIS, 2MASS, ISOCAM and other data. We then compare -/V H ,x and Aj for each 
object, up to unprecedently high extinction. For the p Oph dark cloud with a relatively large sample of 20 bona-fide 
Class III sources, we probe the extinction up to Aj < 14 (Ay < 45), and find a best-fit linear relation Ns,x/A] = 
5.6 (± 0.4) x 10 21 cm" 2 mag" 1 , adopting standard ISM abundances. The other regions reveal a large dispersion 
in the JVh,x/ Aj ratio for each source but for lack of adequate IR data these studies remain limited to moderate 
extinctions (Aj < 1.5 or Ay < 5). For p Oph, the Nn,x/Aj ratio is significantly lower (> 2a) than the galactic 
|' value, derived using the standard extinction curve (Ry = 3.1). This result is consistent with the recent downwards 

1~ \ revision of the metallicity of the Sun and stars in the solar vicinity. We find that the p Oph dense cloud has the 

same metallicity than the local ISM when assuming that the galactic gas-to-dust ratio remains unchanged. The 
difference between galactic and local values of the gas-to-dust ratio can thus be attributed entirely to a difference 
, in metallicity. 



Key words. X-rays: stars - Stars: pre-main sequence - ISM: clouds - ISM: dust, extinction - open clusters and 
associations: p Oph, Cha I, RCrA dark cloud, the Orion nebula, NGC 1333, IC348 

1. Introduction rive M gas from Nr., and some modeling of dust grains 

to go from A v or Aj to M dus t- Although M gajS /M dust a 

One of the key parameters to study galactic evolution, pnQm depends both on the length and direction of the 

star formation, and other fundamental problems in as- line _ of _ sight; the we n_ known remarkable result is that this 

trophysics is the so-called "gas-to-dust ratio". The only ratio hag almogt the same yalue (with gome dispersion ) 

way to_ measure this ratio is to start from the^ "gas-to- throughout the GalaxVj with Mgas /M dust ~ 100. 
extinction" ratio, expressed in terms of the ratio of the 

gas (i.e., hydrogen in all forms) column density N H to- Perhaps the most widely known method to deter- 

ward some sources, to the dust extinction in front of the mine the N H /A V ratio is to use the Lya absorption 

same source, usually expressed in near-IR or optical mag- to measure the H I column density in front of hot, 

nitudes Aj or Ay. The dust extinction is generally derived UV emitting OB stars, and use their spectral types 

from data covering a limited wavelength range in the UV and B - V color excess to derive the extinction. Such 

to near-IR bands. The actual ratio by mass (M gas /M dust ) measurements have been performed in the p ast using 

can then be obtained from the N B /A V ratio, at the cost the Covernicus feohlin. Savage fc Drakeill97Sl) an d WE 

of some knowledge of the physical state of the gas to de- satell ites (|Shull fc van Steenbergj |l985t |Diplas fc Savagj, 

Il994|) . Copernicus also provided a complementary measure 

Send offprint requests to: My Ha Vuong of the H 2 column density from the rotational and vibra- 

e-mail: vuong@discovery.saclay.cea.fr tional levels of hydrogen in the UV. However, the method 
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Table 1. Methods for determining Nji/Ay ratios. 



Method 


A\r max 


Medium 


Length scale 


Ref. 


B stars (UV) 


~2 mag 


diffuse local 


~ few 100 pc 


Bohlin. Savaee. & Drake (1978 ) 


B stars (UV) 


~2 mag 


diffuse local 


~pc 


Whittet 11981} 


SNR (X-rays) 


-~30 mag 


diffuse galactic 


~kpc 


Gorcnstcin (1975) 


SNR (X-rays) 


~30 mag 


diffuse galactic 


~kpc 


Rvter et al. ( 1975) 


SNR + compact objects (X-rays)* 


-~20 mag 


diffuse galactic 


~kpc 


Predehl & Schmitt (1995) 


Diskless T Tauri stars (X-rays)* 


^50 mag 


dense clouds 


~ sub—pc 


This work 



* X-ray spectroscopy with imaging X-ray satellites 



is limited to low extinctions since UV photons can only 
penetrate up to Ay < 2 — 3. 

In this paper, we measure hydrogen column densities 
using X-ray absorption toward young stars. The method is 
conceptually similar to the classical method based on Lya 
absorption in the UV, with two major improvements : 

1. X-ray absorption measures the total hydrogen column 
density (solid state, molecular, atomic and ionized), 
contrary to Lya absorption which measures only the 
atomic hydrogen column density. 

2. X-rays can probe much higher extinctions (with the 
single assumption of a metallicity ratio) than the UV 
range: several tens of magnitudes or more (depending 
on the X-ray energy range available), as opposed to 
only a few magnitudes for the UV, which is ideal for 
the study of dense material. 

More precisely, X-rays are absorbed along 
the line of sigh t by photoelectric e ffec t on 
heavy atoms {Morrison fc; McCammonl . [l 983; 
IWilms. Allen, fc McCravi l200(fl . The X-rav extinc- 
tion cross section varies as E v 2,5 hence rapidl y decreases 
at high X-ray energies (e.g., iRvteii Il996|) . In dense 
regions such as the cores of molecular clouds or the 
immediate vicinity of young stellar objects (YSO), 
self-shielding effect may become important giving a 
reduced contribution to the X-ray absorption. Except 
for this effect which concerns only the relatively minor 
K-shell edges of refractory elements, partially ionized, 
atomic, molecular and solid phases contribute equally 
to the X-ray absorption. In other words, the X-ray 
absorption is dominated not by the H atoms themselves 
but by metals, mainly C, N, and O, and is almost 
insens itive to the presence of grains al ong the line-of-sight 
(e.g., iMorrison fc McCammonl Il983h . 

The method of using X-rays to determine the total hy- 
drogen column density is not new. It has been applied as 
early as 1975 from X-ray spectra obtained usi ng scanning 
X-ray colli mators aboard rocke ts or satellites. iGorensteml 
(1975) and iRvter et al. I l|l975h derived relations between 
interstellar X-ray absorption and optical extinction, from 
the correlation of measurements of supernova remnants 
(SNRs) and the optical E(B — V) color excess, up to 
Ay ~ 30. More recently, IPredehl fc Schmittl l)l995l) used 
ROSAT images to study X-ray halos produced by dust 
scattering around bright compact point sources and super- 
nova remnants, and determined with this technique dust 



extinctions up to Ay < 20. All these studies have assumed 
solar abundances for the ISM absorption. 

However, all the above methods probe mostly the low- 
density, diffuse interstellar medium. High extinctions are 
reached in these X-ray studies only because the line-of- 
sight is long, typically on the kpc scale, so that many 
local effects are averaged out, not only in depth but also 
in projection on the plane of the sky (in the case of SNRs, 
which are extended X-ray sources). 

In the present paper, we want to address this issue for 
the first time in the case of nearby (d < 500 pc), dense 
molecular clouds, where the line-of-sight is dominated by 
H2. The high extinction is here the result not of a long 
line-of-sight, but of the high densities of molecular clouds 
(n > 10 5 cm -3 ); the relevant length scale is here small, 
on the order of parsecs or less. The idea is to make use of 
the fact that "young stellar objects" (YSOs) are ubiqui- 
tous X-ray emitters, and at the same time may suffer high 
extinction, either because they are deeply embedded in, 
or because they are located behind, dense material. We 
show below that after having defined appropriate selec- 
tion criteria, we are able to construct a reliable sample of 
X-ray "candles" allowing the derivation of A^h (from X- 
ray absorption) and Aj (applying standard methods) for 
each object, up to Aj ~ 15 (equivalent to Ay w 45), and 
potentially even more. This is now possible because of the 
spectroscopic capabilities of the new generation of X-ray 
observatories ( Chandra, XMM-Newton) with much higher 
sensitivity and angular resolution in a wide (0.5—10 keV) 
energy band. 

Table ^ summarizes the above mentioned approaches, 
for comparison with that used in the present paper. 

The origin of the X-ray emission from T Tauri stars 
(TTS) is an enhanced solar-type activity generated by 
powerful magnetic reconnect ion flares (see review by 
iFeigelson fc MontmerleL Il999|) . X-ray emission is caused 
by thermal bremsstrahlung and emission lines from an op- 
tically thin coronal plasma confined by magnetic loops, at 
temperatures T x ~ 10 6 - 10 8 K (~ 0.1-10 keV). The av- 
erage Lx/iboi for T Tauri stars is - 10~ 4±1 , or 10 2 - 10 4 
times higher than the active Sun, enabling their detection 
up to large distances. 

Obtaining Ay measurements for the dust component 
of highly obscured young stellar objects can be difficult 
because the y-band flux is often not measurable and 
the infrared (IR) spectrum may be dominated by lo- 
cal circumstellar emission. For this reason, we select a 
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Region name Position (J2000) Obs ID Observation date Exposure time Ref. 



p Oph cores B+F 


16 h 27 m 17 s -24°34'39" 


635 


2000-04-13 


lOOks 


Imanishi ct al. (2001a,b, 2002) 


p Oph core A 


16 h 26 m 34 s -24°23'28" 


637 


2000-05-15 


lOOks 


Daniel. Gaene. & Skinner (2000) 


p Oph mosaic 


16 h 25 m 32 s -24°25'03" 


618 


2000-06-22 


5ks 


This work 


p Oph mosaic 


16 h 28 m 26 s -24°17'33" 


623 


2000-02-08 


5ks 


This work 


Cha I 


ll h 10 m 12 s -76°34'30" 


1867 


2001-07-02 


70ks 


This work 


R CrA 


19 h 01 m 50 s -36°57'30" 


19 


2000-10-07 


20ks 


This work 


IC348 


03 h 44 m 30 s +32°08'00" 


606 


2000-09-21 


50ks 


Preibisch & Zinnecker (2001. 2002) 


NGC 1333 


03 h 29 m 05 s +31°19'19" 


642 


2000-07-12 


50ks 


Getman et al. (2002) 


Orion 


05 h 35 m 15 s -05°23'20" 


18 


1999-10-12 


47ks 


Feieelson et al. (2002) 


Orion 


05 h 35 m 15 s -05°23'20" 


1522 


2000-04-01 


36ks 


Feieelson et al. (2002) 



sample of X-ray emitting stars without circumstellar en- 
velopes or disks. In the contemporary classification scheme 
of pre-main sequence evolution, such stars are identified 
as Class III sources in the near-IR, or ' weak-l i ned' T 
Tauri stars (WTTS) in the optical (e.g., iLadal Il99lj) . 
The WTTS are characterized by the small equivalent 
width of their Ha emission lines (EW(Ha)< lOA) (e.g., 
iMartin et al. I . Il998() . The spectral energy distributions 
(SED) of Class III sources is consistent with a simple red- 
dened blackbody photosphere. 

In this paper, we present a study of the foreground ex- 
tinction of the bona-fide Class III stars in several nearby 
molecular clouds based on Chandra and XMM-Newton ob- 
servations and IR data. We describe the determination of 
the gas column density from the X-ray spectra (0 and 
the dust extinction from the IR data (Jill f° r various re- 
gions. Next, we compare the gas column density and the 
IR extinction leading to A'h.x/Aj (2J, comparing with 
the galactic value and considering the spatial dependence. 

discusses the effects of grain properties, as well as that 
of new stellar and ISM abundances. Conclusions and im- 
plications are presented in ^ 

2. Determination of the hydrogen column density 
from X-ray absorption 

2.1. Observational context 

The principle of deriving the hydrogen column density 
from the X-ray spectrum of any X-ray source is as fol- 
lows. A number of parameters can be adjusted to fit the 
observed spectrum with an input spectrum, which reflects 
the intrinsic X-ray emission mechanism of the source, in 
our case bremsstrahlung continuum plus emission lines 
from an optically thin, hot plasma (10 6 — 10 s K). The 
main fitting parameters are the absorption hydrogen col- 
umn density (-/Vh,x) constrained mainly by the soft energy 
part of the spectrum (below 1—2 keV) as a consequence 
of the energy dependence of the extinction cross section 
discussed above and the temperature Tx of the emitting 
plasma, which is determined by the hard part of the spec- 
trum. It is not unusual for T Tauri stars that the spec- 
tral fit is improved by using two-temperatures plasma to 
describe the continuous temperature distribution of the 



corona. In this paper, we are only interested in deriving 
A^h.x for each source, but the accuracy of this derivation 
depends on the overall quality of the spectral fit, which 
implies a correct modeling of the plasma temperature. 

We consider six nearby star forming regions ob- 
served with Chandra and XMM- Newton: R CrA (~ 
130 pc. iMarraco fc Rvderenl Il98l|) . p Oph (140 pc, see 
Grosso ct al., 2000, and references therein ), Cha I ( 160 pc , 
Whittet et all Il997l). TC 348 (310pc Irlerhiel ll99Sl) . 
NGC 1333 (320 dc. Ide Zeeuw et all Il999l) and the Orion 
Nebula Cluster (450 pc). Part of the p Oph, NGC 1333, 
IC 348 and Orion Chandra observations have been al- 
ready p ublished by Imanish i et al. (2001a) , Get man et alJ 
i|2002l) . IPreibisch & Zinneckerl l(2002ll and iFeigelson et all 
1)20021) . respectively. These papers do quote the gas col- 
umn density Nn but no information on the goodness-of- 
fit is provided, we therefore reprocessed the corresponding 
archived data. 

To compare with previou s X-ray work, we a ssume 
standard solar ab u ndances llAnders fc Ebiharal Il982t 
lOrevesse fc Au"cWsl \m& ICrevesse fc Sa.uvall Il998h as 
used by default in XSPEC for the ISM. We defer the effect 
of adopting more recent abundance determinations to |21 

2.2. Chandra data reduction 

All of the X-ray Chandra observations we use here are 
obtained with the ACIS-I array consisting of four front- 
side illuminated CCDs w ith a field of view about 17' x 17' 
l|Weisskopf et all l2002|) . The spectral resolution E/AE 
is ~ 20. The pixel size is 0'.'5 and the on-axis effective 
area is 600 cm 2 at 1.5 keV. We start the data analy- 
sis with the Level 1 processed events provided by the 
pipeline processing at the Chandr a, X-ray Center. We ap- 
ply the techniques developed by iTownslev et all (l200(t) 
to correct the energy and grade of data events affected 
by the degradation of charge transfer inefficiency (CTI) 
in the energy range 0.5—8 keV. We remove bad pixels, 
events arising from cosmic rays and detector noise using 
the ASCA grades (0,2,3,4,6) filter. We keep events flagged 
as cosmic ray afterglow events for source spectral analysis 
to avoid loosing X-ray events on bright X-ray sources. We 
use the rmf files (response matrix files) and the quantum 
efficiency uniformity (QEU) files provided with the CTI 
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GY135 

N_H,X = 2.0E22 T_X=2 keV L_X = 1.3E30 erq/s 
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Fig. 1. Sample of Chandra X-ray spectra and their corresponding IR SED fitting spectra. The left panels show X-ray 
spectra for three moderately absorbed p Oph sources (GY135, VSSG11 and WL5). Their corresponding IR spectral 
energy distribution (SED) are in the right panels. Note the physical parameters (Nn,x, 2x, Lx) derived from the 
X-ray spectra, using the default solar ISM abundances in XSPEC (§2), and the corresponding (Aj, Tbb, Lhoi) derived 
from the IR SED fitting (§3). See Tables El and H for the errors. 



corrector package 1 . These QEU files are used by mkarf 
in the CIAO package 2 to generate arf files (ancillary re- 
sponse files) consistent with our CTI corrected events. A 
more detailed discussion of these data analysis methods 



1 http:/ /www. astro. psu.edu/users/townsley/cti/install. html 

2 http://cxc.harvard.edu/ciao/ 



is pro vided by iFeigelson et alJ (|2002l) and iGetman et alJ 
ll2002l) . 

We extract events (i.e., individual photons) in the en- 
ergy range 0.5— 8keV for source analysis using circular 
apertures centered on the source position. The events are 
grouped into spectral bins of 15 photons using the ftools 
command grppha. To have enough statistics for the spec- 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 

H-K K-L 



Fig. 2. JHK and JKL color-color diagrams for X-ra y selected (> 100 cou nts) Class III sources in p Oph. The 
solid lines are the intrinsic color s of dwarfs and giants jBessell fc Brettl Il988l) . The dashed lines extending from the 
extremity of the ZAMS/RGB are lRieke fc Lebofskvl ljl985|) reddening vectors representing Aj=5 mag (Ay ~ 16 mag). 
The typical error bar is shown in the upper left corner. 



tral fitting, we keep only Class III sources detected in X- 
rays with a good signal-to- noise i.e., at least 100 counts. 
We measure the background from a neighboring area de- 
void of X-ray sources. Lightcurves are studied for each 
source. For the sources displaying obvious X-ray flares dur- 
ing the observation, we retain for the spectral analysis only 
the quiescent phase, since plasma a bundance changes ar e 
seen to occur in YSO X-ray flares l|Kamata et allll99^) . 
correlated with the rise of temperature. 



The spectra are background subtracted and fitted with 
bremss trahlung, contin uum models using XSPEC version 
11.0.1 l)Arnaudl ll^jlSj)- The fitted models are a combi- 
nation of a n absorption model usin g the atomic cross- 
sections of iMorrison fc McCammonl (|1983) and an op- 
tically thin hot plasma with bremsstrahlung and emis- 
sion lines using the mekal code. The free parameters 
of the fit are the plasma temperature Tx and absorb- 
ing gas column density -/Vh,x- The errors are given for 
90% confidence intervals obtained using the command 
error or steppar in XSPEC. The coronal plasma abun- 
dance is also a parameter of the fitting (see discussion 
in §2.6). Our 1-temperature modeling is relevant if the 
Q— probability, (i.e., the probability that one would ob- 
serve the xti or a larger value, if the assumed model is 
true, and the best-fit model parameters are the true pa- 
rameter values), is greater than 5%. Otherwise we add 
another plasma temperature component to satisfy the cri- 
terion of Q— probability > 5%. 



Fig- HI illustrates a sample of moderately absorbed X- 
ray spectra for three p Oph sources. Table |3] presents the 
spectral parameters and the deduced X-ray properties of 
Class III sources observed with both Chandra and XMM- 
Newton satellites. 



2.3. Chandra observations of the p Ophiuchi dark 
cloud 

The p Ophiuchi dark cloud is a nearby low-mass star 
formation molecular cloud which has been the sub- 
ject of numerous X-ray studies. Early observations with 
the Einstein Observatory satellite showed that YSOs 
are bright and v ariable X-ray emitters i n the 0.2—4.5 
keV e nergy b and llMontmerle et all Il983j) . iKovama et alJ 
l|l994|) and iKamata et alJ (|l997|) detected hard X- 
rays from T Tauri stars w i th A SCA i n the range 
0.5-10 keV. ICasanova et all l|l995l) and ICrosso et all 
l|l997l l2000l) reported deep ROSAT Position Sensitive 
Proportional Counter (PSPC) and the follow-up ROSAT 
High Resolution Imager (HRI) imaging of the p Oph cloud 
dense cores A and F. 

The Chandra observations were made with the ACIS— I 
CCD array in the 0.2-10 keV band. A large area of the p 
Oph region has been covered (~ 40' x 40'). It consists in 
a lOOks exposure observation of cores B and F, a 100 ks 
exposure observation of the core A and seven 5 ks shallow 
exposures forming a mosaic (see TablejUfor d etails) . In the 
core B and F region. Ilmanishi et 

ID l|2001ah identified 58 
X-ray sources with near-IR counterparts. From this data 
set, we select 6 T Tauri stars without m id-IR excess, i.e. , 
witho ut accretion disks (ISOCAM data. IBontemps et all 
2001). We re-analyze the X-ray spectrum of these sources 
to obtain the goodness-of-fit parameters. Five sources have 
enough counts for our study. For the other fields, we iden- 
tified 8 bona-fide Class III sources in the p Oph A field and 
3 Class III sourc e s in th e mosaics using the classification of 
IBontemps et alJ l)200ll) . The selected p Oph sources have 
166 to 7400 X-ray counts. 

Fig. presents the color-color diagrams J — H vs. 
H — K and J—K vs. K — L of our resulting sample of X-ray 
emitting Class III sources in p Ophiuchi with available L 
band photometry. These NIR sources previously classified 
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Fig. 3. JHK and JKL color-color diagrams for X-ray sources in Cha I, NGC 1333, IC348 and Orion Nebula. The 
solid lines are t he unreddened ZAMS an d red giant branch. The dashed lines extending from the extremity of the 
ZAMS/RGB are lRieke fc Lebofskvl l)l985|) reddening vectors representing Aj=5 mag (Ay ~ 16 mag). The plus signs 
(+) are sources representing IR L band excess, open circles (o) are faint X-ray sources (< 100 counts) without IR 
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Table 3. X-ray properties of Class III sources in p Oph, IC348, R CrA, Cha I and Orion Nebula observed with 
Chandra and X-ray properties of Class III sources in p Oph observed with XMM-Newton (using the default standard 
solar ISM abundances in XSPEC). 
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a degrees of freedom. 



b Absorption corrected luminosity Lx (10 30 erg.s 1 ) from the X-ray spectrum in the 0.5—10 keV range. 
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as Class III using ISOCAM data l)Bontemps et alll200l|) 
are located within the area defined by the reddening vec- 
tors. Since ISOCAM data are not available for the other 
regions, we then use similar diagrams to define bona-fide 
Class III sources as those, taking into account their color 
uncertainties, which lie in the area between the redden- 
ing vectors extending from the extremity of the ZAMS in 
both diagrams (cf., Figs|SJ|. 

2.4. Chandra observations of other star forming clouds 

The Chamaeleon complex is one of the star formation 
regions nearest to the Sun (between 160 and 180 pc, 
IWhittet et allfl997j) . It consists of three main dark clouds, 
designated Cha I, II and III. Cha I is the most active and 
hence the most studied dark cloud of the Cha complex. X- 
ray observatio ns of the Cha I region h ave been obtained 
with Einstein llFeieelson fe Krisd. Il989^) and with ROSAT 
l|Feigelson et all Il993) . We use a 70ks exposure archival 
image of Cha I obtained with Chandra (see Table|5J). Eight 
X-ray sources are previously known WTTS and have pub- 
lished JHKL photometry l|Kenvon fe Gomezl EoOlj) . To 
select the TTS without circumstellar disks, we used two 
color-color diagrams (J—H vs. H—K and J—K vs. K—L) 
shown in Fig. [3] and apply the same criterion as for p Oph 
(Fig.EJ). We then select 4 X-ray emitting Class III sources 
with more than 100 counts. 

Like Cha I, the R Coronae Australis (R CrA 
cloud) has a lso been extensively studied in X-rays : 
with Einstein Walter fc Kuhil . ll98lUWalter et allll997h . 
with ASCA faoyama et all Il996h and with ROSAT 
ijNeuhauser et alll2000(l . In the 20 ks archival Chandra ex- 
posure, three Class III sources are identified but only two 
are sufficiently bright for our study. 

NGC 1333, an active star forming region within the 
Perseus molecular co mplex, wa s prev iously observed and 
st udied in X-rays bv IPreibischl l|l997f) using ROSAT and 
bv lCetman et al.1 l 2002h with Chandra. In the 37.8 ks ex- 
posure time, they detected 127 X-ray sources and iden- 
tified 95 with previously known cluster members in the 
optical and near IR. Spectroscopy is not available for 
these known cluster members. Seventeen X-ray sources are 
ident i fied with JHKL sources l)Aspin. Sandell. fc Russell 
119941 lAsnin fc Sauddl ll 997ft and have been plotted in 
the two color-color diagrams (Fig. EJ. We find 14 sources 
within the area defined by the reddening vectors extend- 
ing from the extremity of the ZAMS in both diagrams in- 
cluding the foreground B star BD+30°547. These sources 
are symbolized as circles in both diagrams. Three sources 
meet the X-ray selection criterion, but unfortunately they 
fall outside of thi s area, so that their st atus is unclear 
(see discussion in lAspin fc Sandell Il997|) . We have not 
included these sources in our sample, and therefore the 
whole region, in spite of its interest, has to be excluded 
from our study for lack of adequate data. 

IC 348 is a nearby star forming region also be- 
longing to the Perseus molecular cloud complex. 



The most massive cluster member is the B5 V star 
BD+31°643. ROSAT X-ray ob servat ions of IC348 
l|Preibisch. Zinnecker. fc Herbid . Il99fih had shown 
that this region co n tains over a hundred stars. 
IPreibisch fc Zinneckerl l|200ll l2002h obtained a 50 ks 
exposure image of IC348 with ACIS— I. They detected 
215 X-ray sources and identified 115 of these sources with 
known cluster members. They used the absorption de- 
rived from optical/near- IR spectroscopy I Luhman et all 
ll99Sh to fit the spectra instead of deriving it from 
the X-ray spectra. Therefore, we have re-analyzed this 
dataset to derive the gas column density from our X-ray 
spectra. Op tical /near-IR spectro s copy of this region 
is available llLuhman et all ll 9981 lLuhma.nl ll 999ft . But 
IPreibisch fc Zinneckerl l|2002h noticed that many TTS 
selected on the basis of Ha equivalent width EW(Ha) < 
10 A have a JHK infrared excess. We first pre-select the 
X-ray sources whose counterparts have EW(Hq) < 10 A 
and we found 29 sources. N ext, we use JHK phot o metry 
and the 3.5 /im data from lHaisch. Lada. fc Ladal l|200ll) 
of these sources to construct the color-color diagrams 
(Fig. Twenty X-ray sources are located within the 
reddening band in both color-color diagrams which have 
no L band excess. Eight have > 100 counts required for a 
good fit to the X-ray spectrum. 



The Orion Nebula is the most distant highly active 
star forming region of our sample. Observations with the 
Einstein, ROSAT and ASCA satellites have been car- 
ried out but could identify only a small fraction of X-ray 
sources of the cluster population due to sensitivity, angular 
resolution and bandwidth limitations. With much higher 
sensitivity and angular resolution, Chandra recently de- 
tected more than 1,000 sources in the ONC identifie d 
with previously kn own members ijFeieelson et all [2002) . 
iHillenhraudl l|l997l^ studied the optical spectroscopy of 
~ 900 stars in this region but did not provide informa- 
tion on the Ha equivalent width (to distinguish between 
CTTS and WTTS), due to the contamination of this 
line by the emission of the blister H II region associated 
with the Orion Nebula. Using JHKLMNQ catalog in 
iHillenbrand et al.1 l)l998l) . we corr elate this catalog with 
the X-ray source list detected by iFeigelson et al.l (|2002) 
and find 16 sources which have JHK and L data and 
a good X-ray spectrum. We select seven of these sources 
which have no L-band excess and are located within the 
reddening band (Fig. |HJ). Six s ources have small or zer o 
if-band excess A(I-K) < 0.3 Hillenbrand et allll998h 
The c olumn density derived from the IFeigelson et all 
(2002^) paper is not sufficiently accurate because of the 
automated fitting procedure of a large number of sources. 
We re- analyze the X-ray spectra of these 6 sources and find 
improved spectral parameters which are listed in Table 01 
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Fig. 4. Confidence contours of JVh,x versus temperature spectral fit of IRS55 in the p Oph cloud detected by both 
Chandra (left) and XMM-Newton (right). The contours are la, 2a and 3<r levels respectively. The values of Nu,x from 
these two independent observations are consistent within la confidence. 



2.5. XMM-Newton observation of the p Ophiuchi 
dark cloud 

The p Ophiuchi cloud was also observed with XMM- 
Newton in the 0.2-10.0 keV band for 32 ks on 
2001 February 19 centered at R.A. = 16 h 27 m 26.0 s 
and DEC=-24°40'48.0". XMM-Newton features three 
coaligned X-ray telescopes with 6" FWHM angular reso- 
lution at 1.5 keV an d with large total eff ective area up to 
5000 cm 2 at 1 keV l|jansen et all l200l|) . The European 
Photon Imaging Camera (EPIC) consists of two MOS 
(MOS 1 and MOS 2) and one backlighted PN, X-ray CCD 
arrays placed on each focal plane of the telescopes. The 
EPIC provides capability of imaging with a field of view 
of 30' and spectroscopy with moderate spectral resolution 
of - 60 eV at 1 keV. 

The data reduction followed standard procedures using 
the XMM-Newton Science Analysis Software (SAS version 
5.3). We ran emchain and epchain to obtain photon event 
lists. We removed bad pixels and selected single, double, 
triple, and quadruple pixel events for MOS and single and 
double pixel events for PN using evselect. We extracted 
light curves from whole field of view and excluded the time 
region where the X-ray counts are extremely high, which 
could result from solar flare protons. 

The background-subtracted X-ray spectra were ob- 
tained from 8 Class III sources with good statistics (150- 
3200 counts). The counts were extracted from the 90% 
enclosed energy circle. In the case other X-ray sources are 
found close to the Class III sources, we adjusted the radius 
of the reference regions to minimize the contamination 
by these sources. The background X-ray counts were ex- 
tracted from annular regions around each Class III source. 
The X-ray spectra were obtained from each EPIC CCD 
detector (MOS1, MOS2, and PN), and fit simultaneously 
with the same model used for Chandra spectral analysis 
using XSPEC. 



We find that three sources (IRS 55, SR12 and WL 19) 
of our sample are in common with Chandra. This allows a 
test for possible systematic errors between the two satel- 
lites. 

2.6. Reliability of the X-ray derived iVn.x values 

We consider here several contributions to the uncertainty 
of the -/Vh,x values derived from the Class III X-ray spec- 
tra: possible systematic errors associated with satellite cal- 
ibrations; statistical errors from the limited source counts; 
possible systematic errors associated with the uncertain 
X-ray spectrum of the underlying star; and with uncer- 
tainties in the stellar coronal abundances. 

The spectra are fitted using x 2_u t statistics with 
the minimization technique of Levenberg-Marquardt 
l|Bevington fc Robinsori 119921. The errors of 7Vh,x from 
the model fitting are given for 90% confidence intervals 
obtained using the command error in XSPEC. Table [21 
reports these values. 

We first compare the instrumental errors of the de- 
termination of iVn.x between Chandra and XMM-Newton 
for three Class III sources observed by both satellites. 
Fig. 0| shows the confidence contours of Ah,x and two- 
temperature fits for IRS55 in p Oph. The agreement be- 
tween the values of Ah,x from these two independent ob- 
servations, which are consistent to better than 1 a, is re- 
markably good, and gives us confidence in the reliability 
of the X-ray spectral fitting procedure. 

Another factor of uncertainty may be the influence of 
statistics, i.e., whether t he number of c ounts affects the 
determination of Ah,x- iGetman et al.l l|2002|) have per- 
formed simulations with sources of known spectra us- 
ing the f akeit utility in the XSPEC package. Since our 
sources are physically of the same nature (T Tauri stars), 
their results are also valid in our case. They found that 
the standard deviations of the derived column density val- 
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Fig. 5. Comparison of the v alues of the dust extinctio n A] 
derived from spectroscopy l|Luhman k. Riekel Il999|) and 
those derived from the IR SED fitting (this work) for 10 
Class III sources in p Oph. The dashed line is the linear 
fitting line. The solid line is the equality line A} (this work) 
= A] l|Luhman fc Riekel. 1 1999j) . 



ues AiVjj.x (x 10 22 ) range from ~0.6 for sources with < 
300 counts to ~0.3 for sources with > 300 counts. These 
values are included in the error bars given in Table |3| 

We have also checked the effect of the coronal abun- 
dance on the value of iVn,x- The coronal plasma abun- 
dance is first fixed at the best fit value. We then test with 
different values for each spectrum until the % 2 deteriorates 
beyond a Q— probability = 5%. We find that changing the 
coronal abundance by a factor of ~ 2 leads to a < 10% 
variation in the value of iVj^x, a range which is already 
comprised in the error bars. 

Lastly, a systematic uncertainty in ATh.x could also 
arise due to a possible unobserved, absorbed soft compo- 
nent. Such a component would be a priori undetectable 
in case of high extinction, but could modify the spectral 
fit and surreptitiously affect the determination of -/Vh,x- 
We have checked this by adding a fake soft component to 
the X-ray spectra. We have used the observed spectrum 
of WL5, which has the highest iVn,x and enough counts 
to study the effect of the soft component. Using fake it 
utility in the XSPEC package, we simulated the spectrum 
of WL 5 by adding a soft component at 0.3 keV to a 
hard component at 2 keV. We fitted the resulting spec- 
trum with one-temperature model as seen with the real 
spectrum. We found that if there is a soft component in 
the spectrum of WL5 brighter than 5% of the total flux, 
it would be picked up by XSPEC. This is already what 
we find in other sources which require a two-temperature 
model (like IRS 55 just mentioned, see Fig. @J|. For WL5 
we conclude that no significant unseen soft component can 
affect iVn,X) an d a fortiori the other sources with smaller 
Nh,x are not affected. 



Table 4. 


Comparison 


between 


A j , and our 


SED 


blackbody 


temperature 


fitting 


with those 


from 


Luhman & Riekel ( 19991) for the p Onh sources. 






SED fitting^ 


K — band spectroscopy^' 


Sources 


Aj 


Tbb 


4j 


T e fl 


GY12 


5.02-5.23 


3549-3925 


4.6 


3955 


GY135 


3.61-3.64 


4323-4690 


4.4 


4205 


GY156 


5.05-5.22 


4268-4733 


6.1 


3955 


LFAM8 


5.51-5.81 


3434-3873 


6.6 


3595 


VSSG11 


3.50-3.51 


3894-4262 


4.0 


3850 


B162607 


4.98-5.11 


3533-4103 






SI 


3.57-3.81 17665-19935 




18700 


DoAr21 


1.80-1.97 


4614-4929 


1.6 


5080 


SR20 


2.47-2.53 


3677-3945 


2.1 


3955 


GY289 


6.39-6.81 


3910-4420 






IRS55 


1.41-1.51 


3950-4175 






WL5 


12.26-13.93 


5775-6970 






SR12 


0.21-0.22 


3610-3640 


0.20 


3850 


GY112 


0.96-1.05 


3505-3745 






GY193 


1.80-1.87 


3360-3650 






GY194 


2.19-2.24 


3550-3865 






GY253 


7.08-7.59 


3350-3820 


8.2 


3350 


GY410 


2.51-2.53 


3870-4120 


2.7 


4060 


B162730 


2.44-2.49 


3555-3850 






WL19 


14.46-16.45 


6000-7210 







(a) This work; (b) Luhman & Riekc (1999) 



3. Determination of the extinction from IR 
photometry 

3.1. IRSED fitting 

Different methods may be used to derive the extinc- 
tion by dust grains. Because the spectra of Class III 
sources lie in the near-IR range, the extinction Aj is di- 
rectly derived, instead of the more usual Ay. Available 
optical/IR spectroscopy leads to the spectral type de- 
termination and intrinsic colors which, when compared 
to the observed colors and an extinction law, give an 
estimate of the dust e xtinction. This is the case for 
most sources in IC 348 llLuhman et all 1998; iLuhmanl 
199 9^ Cha I jLawson. Feigelson. fc Huenemoerderl 11996: 
lOauvin fc StromL Jl992lL R, CrA llWalter et all Il997lb 
QNC |li ill^nbrancijl997j) and for some objects in p Oph 
(Lu hman fc Riekel [1999J). These sources generally have 
low or moderate extinction, allowing optical/IR spec- 
troscopy. 

If the spectral type is unknown, the dust extinction 
is determined by fitting the near to mid-IR spectral en- 
ergy distribution (SED) with a blackbody spectrum. This 
is the case of some embedded sources in p Oph where 
Aj can be as high as 14 magnitudes (Ay — 50 mag) 
and optical/near-IR spectroscopy cannot lead to an ac- 
curate determination of the spectral type. We use the / 
band {DENIS), JHK photometry (2M ASS) and 6.7 urn 
and 14.3 /im photometry from ISOCAM l|Bontemps et all 
to construct the SED. For some sources, we also 
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have £(3.5/iml M (4.8/im) and N(10fim) photometry 
^Greene et all Il994^ . 

The wavelength dependence of the interstellar extinc- 
tion and polarization provides constraints on the charac- 
teristics of interstellar grains. It is well-established that 
the extinction law depends only on the ratio of to- 
tal to selective extinction Ry = Ay /E(B — V) ~ 3.1 
for the diffuse interstellar medium llRieke fc Lebofskvl 
Il985l ICardelli. Clavton. fe Mathisl Il989|) and is higher in 
some regions. An important point is that, as shown by 
ICardelh. Clavton. fc Mathisl l)l989|) . the determination of 
Aj is independent of Ry, contrary to the determination 
of Ay (see below). 

iDrainel l|l989j) reviewed the extinction by interstel- 
lar dust at IR wavelengths and showed that for wave- 
lengths between 0.7 and 7 /im the observed extinc- 
tion, both in the Galaxy and in the Magellanic clouds, 
is consistent with a power-law A\ ~ A~° w ith in- 
dex a ~ 1.75. ICardelli. Clavton. fc Mathisl (|l989j) adopt 
the iRieke fc Lebofskvl lll985h curve with a ~ 1.6. 
iMartin fc Whittetl l|l990l) established a power law extinc- 
tion curve with a ~ 1.8 extending from 0.8 to 4.8 p,m for 
both diffuse and dense clouds. In this work, we adopt a 
= 1.8 but consider the effect of values ranging from 1.6 to 
2.0. Beyond 4.8 p , we use the extinction law given by 
iLutz et all l)l996fl based on the ISO/SWS observation of 
the galactic center. 

We fit the SE D of each object using the absorbed 
blackbody model ijMvers et all Il987|l The photospheric 
flux of the source -Fbb(Ibb) is assumed to be a blackbody 
at temperature Tbb. 



Fobs — F D b(Tbb) x e 



with 



Aj ,A x , 

x (— — 

1.086 y Ay 



A 



j 



1.086 



1.25 ' 



(valid for A > l^m), where Fobs and Fbb are the observed 
and blackbody flux in Jansky, t\ is the optical depth at 
the wavelength A and a is the power law index. We assume 
that the dust emission is negligible. The best estimate of 
Aj is found by using a least-squares fit to the blackbody 
model. The fitting parameters are then Tbb, Aj. We vary 
the power law index a from 1.6 to 2.0 to obtain the re- 
sulting uncertainties on the parameters. 

Eleven of 20 X-ray emitting Class III sources in p 
Oph have IR spectra an d therefore known spectral types 
i Luhman fc Ri ek3. ll999h . These authors derived the ab- 
sorption Aj using near IR photometry and spectral type. 
The nine other Class III sources of p Oph lack IR spec- 
troscopy. We derive a blackbody fit of the SED of these 
objects to estimate the extinction value, which is given in 
Tabled Fig.r^illustrates a sample of three fits correspond- 
ing to the three X-ray sources analyzed in the previous 
section. 



Table 5. The N-^/Ay ratios from the literature. (See also 
Table 1.) 



Ref. 



1.8xl0 21 



(a) 



(a) 
2 (a,b) 
2.0xl0 21 cm" 2 (b) 
1.8xl0 21 cm" 2 (b) 



2.2XKT 1 cm 

l21 



Bohlin, Sa vage, fc Drake (1978) 

Whittet (1981) 
Gorenstein ( 19751 Rvter (19961 
Rvter et al. (19751 
Predehl fc Schmitt (1995) 



(a) UV; (b) X-rays. 



3.2. Comparison with other extinction determinations 

Table 0] compares the parameters (Aj, Tbb) derived from 
the S ED fitting with the values from ILuhman fc Riekel 
l|l999|) . The two values of the extinction and the temper- 
ature represent the range resulting fro m the variation of 
the po wer-law index a from 1.6 to 2.0. lLuhman fc Riekel 
dl999|) used if-band spectroscopy in p Oph to obtain de- 
termine spectral types and extinctions. Our fitting values 
are generally in agreement with those given in the litera- 
ture when available. 

Fig.[3conipares the values of t he dust extinction Aj de- 
rived from iT-band spectroscopy l|Luhman fc Riekel Il999) 
and those derived from the IR S EP fitting (This work). A 
linear relation can be found: Aj [ l|Luhman fc Riekel Il99^ ] 
=1.16 x Aj [This work] - 1.06. The values of Aj from 
spectroscopy is ~ 10% larger than our SED fitting values. 



4. The ./Vh,x/Ai ratio in nearby star forming 
regions 

4.1. Comparison with the galactic value 

Fig. shows the correlation between the hydrogen col- 
umn density derived from the X-ray spectra A^h.Xj and the 
dust extinction Aj derived from optical/IR spectroscopy 
for 19 bona-fide Class III TTS in the p Oph cloud, as- 
suming the default solar ISM abundances in XSPEC. The 
extinction Aj is probed up to 14 mag (Ay ~ 45 mag). 
For the sources where Aj is determined from the black- 
body model, the horizontal error bars represent the range 
in Aj obtained by varying the index a from 1.6 to 2.0 
( t)3.1fl . The resulting best-fit linear relation is 7Vh,x=5.57 
(± 0.35) xlO 21 Aj (cm" 2 mag" 1 ). 

However, the gas-to-dust ratio is usually studied 
with the help of Nn/Ay. In our Galaxy, this ratio is 
(7V H /Av) g ai = 18 — 2.2 x 10 21 cm -2 per visual magni- 
tude (Table EJ). This galactic relation, obtai ned using X- 
ray absorption, and valid up to Ay ~ 30 mag dGorensterri , 



19751 iRvter et al. T ll975tlRvterlll996irpredehl fc Schmittl. 



Il995 , see Table 1) is consistent with the one determined 
using Lya absorption in the UV limited to low extinctions 
(Ay < 2-3 mag). 

We want now to compare these values with that ob- 
tained above from the Nu/Aj correlation. This requires 
the conversion of Ay into A j , which depends on the shape 
of the dust extinction curve. ICardelli. Clavton. fc Mathisl 
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Fig. 6. X-ray absorption column density -ZVh.x versus visual extinction Aj for Class III sources in p Oph (a), RCrA 
and Cha I (b), IC348 (c), Orion Nebula (d) using the default standard solar ISM abundances in XSPEC. The error 
bars in X-ray column densities are represented by vertical bars at 90% confidence (1.6a). The dashed line in p Oph 
is the fitting curve: 2Vh=5.57(± 0.35) x 10 21 Aj (cm -2 ). The shaded area is the range of published galactic gas/dust 
relationship (see Table |3J) ■ The symbol size is proportional to the number of photons detected in X-ray spectra in log 
(from 195 to 7378). Note that Aj is up to 14 mag (Ay ~ 45 mag) in the p Oph region whereas Aj < 1.5 mag (Ay < 
5 mag) in the other regions. 



( 1989J) have shown that this shape, for a variety of in- 
terstellar regions, is well represented with a single pa- 
rameter Ry = Ay/(A-g, — Ay) where A-q is the dust 
extinction in the B-band. From this parameterization, 
they found that A } /Ay = 0.4008 - 0.1187(i? v /3.1) _1 . 
We use this empirical relationship, and the average galac- 
tic value of Ry = 3.1, to convert Ay into Aj, yield- 
ing the galactic gas-to-dust correlation in the J— band: 
(N H /Aj) gal = 6.4 - 7.8 x 10 21 cm" 2 mag" 1 . This range 
is indicated by the shaded area in Fig. and compared 
with our results for p Oph: it is clear from Fig. that 
the Nh/Aj correlation lies significantly below (>2o~) the 
galactic values. 

The four other regions (Cha I, RCrA, IC348 and 
Orion) are probed only to low extinctions, A] < 1.5 mag 
(Ay < 4). Fig. compares the hydrogen column den- 
sity and the dust extinction toward 4 bona-fide Class III 
sources in the Cha I dark cloud. Three points deviate from 



the galactic relation. Fig. and HJ1 represent 8 lines of 
sight toward IC 348 and 5 lines of sight toward the Orion 
Nebula Cluster (ONC). Note that Fig. 01 displays directly 
Ay, instead of Ay. the values of Ay for Cl ass III sources 
in ONC are taken from iHillenbrandl l(l997t ) using optical 
spectroscopy. Due to the small sample in each of these re- 
gions, we cannot obtain a relationship between -/Vh,x and 
A} like in p Oph. We note that t he comparison o f TVh.x 
and Ay has alre ady been done bv iFeigelson et all l|2002|) 
for the ONC and lCetman et all l|2002h for NGC 1333. but 
only in global terms, for instance without distinguishing 
between Class III sources and Class II sources (TTS with 
disks), or with automated fits. In these studies, the gen- 
eral behavior shows a large dispersion, and no correlation 
between JVh and Aj or Ay can be found. For comparison, 
there are only two data points in p Oph having Aj < 1.5, 
so that spanning a large range in A] is crucial to look for 
a correlation. 
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In conclusion, for p Oph we find that (i) a tight correla- 
tion exists between Ah and Aj up to Aj ~ 15, and (ii) this 
correlation lies significantly below the galactic correlation 
obtained by converting Ay into Aj using the canonical 
value Ry = 3.1 and Nn/Ay determined from UV and 
X-ray measurements. For the four other regions studied 
(Cha I, R CrA, IC 348, and the ONC), the data probe 
a too small extinction range (Aj < 1.5), and no similar 
correlation can be found. 

4.2. Spatial dependence of the N-h/Aj ratio 
4.2.1. High extinction: p Oph 

We now investigate whether the deviation from the galac- 
tic iVu vs. Aj correlation found in p Oph, and the dis- 
persion of the TVh vs. Aj or Ay data points, depend on 
the location of the sources within the clouds. Fig. sym- 
bolizes Ah,x/Aj along individual lines of sight in relation 
with optical, mid— IR and molecular maps. Lines-of-sight 
with values close (within 10%) to the Ah vs. Aj correla- 
tion are indicated by circles; those for which the Ah/Aj 
ratio deviates by more than + or — 10% are indicated 
by triangles, pointing upwards or downwards respectively. 
We find no spatial correlation between the value of the ob- 
served -/Vh,x/Aj ratio and the DCO + contours indicating 
the dense cores. 

We now discuss more precisely the two sources with 
the highest Aj , particularly WL 5 which plays an impor- 
tant role in the Ah,x _ Aj correlation: WL5 (Aj ~ 14) and 
WL19 (Aj ~ 16). 

• WL 5: This source is located in the core B and 
has b een extensively studied at many wave lengths: ra- 
dio bv lAndre . Montmerle. fe Feigelsonl l)l987|) . mm obser- 
vat ion bv I Andre'^'^oTrtinerl^~Il994T) . near and mid-IR 
by iGreene et all l|l994f) . IGreene fc Meverl l|l995j) deter- 
mined an F7 spectral type, bolometric luminosity 33.5 
Lq and extinction Ay— 51 mag, using near-IR spec- 
troscopic obse rvations. This value of Ay is overesti- 
mated because IGreene fc Meverl l|l995j) used the conver- 
sion Ail Ay = 0.265, which corre sponds to Ry — 2.7 
i Cardelli. Clavton. fc MathisL Il989l) . The 100 ks Chandra 
observation provided 1322 photons for this source. WL5 
is also detected in the XMM-Newton observation but its 
spectrum is contaminated by a Class II source (IRS 37) 
located only 7" away. Given the large number counts in 
the Chandra spectrum, there is no particular problem, and 
.ZVh.x is accurately derived (Fig.^). This source is included 
in Fig.EJi. 

■ WL 19: This object was class ified as a Class I YSO 
by Wilking. Lada. & Young (1989) and as a Class I I by 
lAndre fc Montmerle! l|l994h . iBontemos et all l|200lj) de- 
tected no mid-IR excess from this source and concluded 
that it is likely a luminous Class III source located be- 
hind the cloud (i.e. not embedded within the cloud) with 
very high extinction. This source was observed by both 
Chandra and XMM-Newton. The 100 ks Chandra spec- 
trum provides only about 120 photon counts and the 30 ks 



XMM spectrum has 150 counts. The spectra are shown in 
Fig. El The values of Ah,x derived from these two spectra 
with low statistics are ~ 2 a apart. For this reason, we 
have not included this source in the Ah,x/Aj relation. 

4.2.2. Low extinctions 

Unlike p Oph, the other star forming regions have only one 
Chandra observation and their spatial coverage is limited 
by the ACIS ficld-of-vicw (17' x 17'). We could build only a 
small sample of Class III sources in Cha I, R CrA, IC 348, 
and the ONC. In these regions, Fig.Eb, c, d show large de- 
viations of the ./Vh,x/Aj ratio (> ±10%) from the galactic 
value. Fig. shows the spatial distribution of our sample 
of Class III sources in IC 348. The background image is 
the DSS / band image, and the contours represent IRAS 
60 pm isophotes, taken as column density tracers. As in 
p Oph, we use circles and triangles pointing upwards and 
downwards to symbolize deviations of the Ah/Aj ratio, 
here with respect to the galactic value. Again, no corre- 
lation with column density is apparent. The reasons from 
the deviations observed in Fig. are unclear, although 
they would probably play a minor role if the available 
range in Aj was larger as in p Oph. 

5. Constraints on the gas and dust properties 
from the JVh,xMj ratio 

The properties of the gas and dust in p Oph can be con- 
strained from our measurement of the gas-to-dust ratio. 
Indeed, our value of the slope of the Ah.x — Aj correla- 
tion in p Oph is lower than in the Galaxy. This may arise 
(i) because we underestimate the total gas column density 
Ah,x (through the metallicity Z) or, (ii) because the grain 
properties (abundance with respect to the gas, extinction 
per unit mass in the J— band) are modified toward the p 
Oph cloud. 

5.1. Effect of metallicity 

As recalled in §I\ the total hydrogen column densities are 
inferred from the photoelectric cut-off observed in the X- 
rays. Such a cut-off depends on the column density of 
metals Az, which is then converted into hydrogen column 
density Ah by assuming a set of element abundances, ir- 
respective of whether the heavy elements, mostly C, N, O 
(which we will refer to as "metallicity" ) are in the gas or 
in the grains. 

Solar abundances (determined from analysis of 
the solar photosphere or meteorites) have gener- 
ally been used as the reference abundance for the 
TSM JAnders fc FhiharaL ITflsj iGrevesse fc And^sl Ffls3 
iGrevesse fc Sauvai Il998j) . and this is what we have used 
so far in this paper: as noted previously, these abun- 
dances are the default in XSPEC. We now refe r to these 
abundances as "old" solar abundances. Recentlv.lHo lwegerl 
l)200ll) : lAllende Prieto. Lambert, fc Asplundl l|200l[ l2002) 
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Fig. 7. Spatial distributio n of our sample of Cla ss III sources in p Oph. Background: DSS R optic al image. Foreground : 
ISOCAM7+15 a m image llAbereel et al.lll996 L This combination background map is taken from lGrosso et alJ lfeoOO'). 



The DCO + map of lLoren. Wootten. fc Wilkind fl990) is shown by contours and the dense cores are labeled. The circles 
(o) are sources with N-r^x/Aj normalized by the regression line (5.57 x 10 21 ) with a correlation ± 10%. The triangles 
up (A) are sources with normalized values of A^x/^j higher than 1.1 and the triangles down (y) are lower than 0.9. 
The dashed lines are Chandra and XMM field-of- views. 



data and folded mode 



data and folded 



WL19nf_bin20.pi 



s27_WU9_pn_bin.pl 
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channel energy (keV) 
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Fig. 8. Chandra (left) and XMM-Newton (right) spectra of WL19 (A H ,x = 7 - 14 x 10 22 cm" 2 , T x 
Lx = 0.4 — 1.5 x 10 30 erg.s -1 , using the default standard solar ISM abundances in XSPEC). 



1.5 - 3 keV, 



published a "new" determination of photospheric so- 
lar abundances, in particular for carbon and oxygen. 
Therefore, we have included these "new" solar abundances 
in XSPEC and re-calculated the derived Ah,x from our X- 
ray spectral fits. We find that for a given Aj, A'h.x is — 
20% higher than when using the "old" solar abundances. 
This is illustrated in Fig. llOb . where we compare the X-ray 
derived hydrogen column densities using the "old" solar 



and the "new" solar abundances. We find: Ah (new so- 
lar) = 1.32 x N u (old solar) - 0.15 x 10 22 cm" 2 . The 
linear relation between -/Vh.x and Aj for p Oph becomes: 
7V H (new solar) = 7.24(± 0.46) x 10 21 Aj (cm~ 2 ). This new 
relation is then in remarkable agreement with the relation 
(A r H Mj) g ai = 6.4 - 7.8 x 10 21 cm -2 mag" 1 quoted above 
for the Galaxy when using the galactic value i?y = 3.1. 
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Fig. 1 0. Comparison between the d e rived -ZVh,x using "new" solar abundances ilHolwegerL l200ll 
lAllende Prieto. Lambert, fc Asplundl . l200lL 1 2002) an d (a) those using "old" solar abondances, (b) those using 
photosp heric F fc G star abundances ([Sofia fc MeverL l200ll iR.eddv et all l2003h . (c) those us ing diffuse ISM abun- 
dances l|Wilms. Allen, fc McCravl l2f)0(t and (d) those using diffuse ISM B-star abundances l|Sofia fc MeverL 120011) . 
The solid line is the equality line. 



Abundance measurements have also been made in the 
local diffuse ISM with the HST toward field B stars 
llWilms. Allen fc McCravl 1200(1 ISofia fc Meverl 1200 ll) . 
More recently. iReddv et all l|2003l) presented photospheric 
abundances of 181 F & G disk stars. We shall refer 
to these abundances as the "diffuse ISM Wilms et al. 
2000", "diffuse ISM B star" and "photospheric F & 
G star" abundances, respectively. These total gas plus 
dust abundances are 20-30% lower than t he "old" solar 
abundances (d ependin g on the elem ent ) (S now fc Witrl. 
199fiUCardelli et allll99rilMever. Cardelli. fc Sofialll997t 
Sofia fc MeverU200lh . Therefore, we have re-run our X-ray 
spectral fits for the p Oph sources using these abundances. 
Figs.llOb. c, d show that the derived X-ray column densi- 
ties using these diffuse ISM and photospheric abundances 
are almost equal (± ~ 10%). Noting that these recent 
measurements probe relatively short lines of sight (< 1 
kpc for B stars and < 150 pc for F & G stars). We shall 
define this new set of abundances simply as the "local" 
abundances. 



We conclude that the lower value of -/Vh,xA4j for p Oph 
may be fully accounted for by adopting the metallicity 
corresponding to the "local" abundances. 

5.2. Effect of grain properties 

However, recall tha t several works 

(e.g., IVrba. Covne. fc TapiaL Il998h showed that in 
p Oph, Ry is not galactic. In this section, we keep the 
"old" solar abundances and discuss the effect of grain 
properties in relation with Ry. 

As shown in Appendix A, the A^.x/^j ratio is pro- 
portional to the gas-to-dust mass ratio M gas /Md us t (here- 
inafter G/D). Thus, the low A^.x/^-j ratio could be ex- 
plained if G/D = 70 (instead of the average galactic value 
of 100). Assuming that t his higher dust mass is due to the 
formation of ice mantles l|Tanaka et al.. 1990), about 80 % 
of the gas phase oxygen should be depleted onto grains 
(with [0/H] gas = 3 1Q- 4 . ISavage fc Semhachl <|l99filn . We 
do not consider this possibility in greater detail however 
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Fig. 9. DSS / band image and spatial distribution of 
our sample Class III sources in IC348 with IRAS 60 [xm 
contours overlaid. The circles (o) represent sources with 
N-a t x/Aj normalized by the galactic value of Nn/Aj with 
a correlation ± 10%. The triangles up (A) are sources 
with normalized values of -/Vh,x/Aj higher than 1.1 and 
the triangles down (v) are lower than 0.9. 



because the depletion in p Oph is not expected to be strong 
(see Appendix B). 

Another possibility to lower the -/Vh,x/Aj ratio is that 
the Aj /Ay ratio is larger. This ratio and the shape of the 
extinction curve for A < 1 /iin depend on Ry. This pa- 
rameter is closely related to the mean size of dust grains 
as shown in Fig. A.l. Changes of Ry thus trace pro- 
cesses affecting the dust size distribution. In particular, 
. Ry is expected to increase in cold, d ense molecular clouds 
I Cardelli. Clayton, fc Mathi a. 119891) where grain coagula- 
tion probably occurs (|Stepnik et all 120031) . 

From polarime t ry in optical bands, 

IVrba. Covne. fc Tapial l|l993h found Ry = 4 in a 
peripheral region of the p Oph cloud where Ay < 
3 mag. To extend this study to regions with larger 
Ay's requires near-IR or mid-IR polarimetry, which 
has never been performed before. To derive Ry in p 
Oph, we assume that Nn^/Ay is constant throughout 



the Galaxy, equal to 1 .9 xlO 21 cm 2 mag' 1 (±20%, 
e.g., iKim fc Martini Il996|) . We thus find Ry = 6.0 ± 2.5, 



value brackets that of lVrba. Covne. fc Tapial l|l993|) . Note 
that because of the functional dependence between Ry 
and Aj/Ay, the uncertainty on Ry is large. Therefore, 
with the present measurement of -/Vh,x/Aj, it is difficult 
to reach a firm conclusion on the value of Ry and on its 
possible variation from the periphery to the center of the 
cloud. 



5.3. Constraints on the dust size distribution in p Oph 

We now use a physical model of gas and dust properties to 
constrain the grain size distribution and the gas-to-dust 
mass ratio in p Oph. The details of the grain modeling arc 
described in Appendix A. 

Using the observationnally derived value of Ry for p 
Oph, including the uncertainties (see H5.2(l . we find the 
minimum grain size a m ; n = 0.025 /im for Ry — 4 and 
0.07 /im for Ry = 6 (see Fig. lA.lf) . These values are ~ 
5—10 times larger than in the galactic diffuse medium. 
Fig. lA.ll shows that the average grain size < a > increases 
with a m j n . For p Oph, it is ~ 0.035—0.095 /im instead of 
8xl0~ 3 /im (Eq. (A. 4)) in the galactic diffuse medium. 
We also find that the corresponding gas-to-dust mass ra- 
tio G/D for p Oph lies between 80 and 95 using the MRN 
size distribution as described in App e ndix A. This is to 



i.e., A j I Ay = 0.34+£^ (instead of 0.28 for Ry = 3.1) 
usin g the parameterization of the ext inction curve given 
bv lCardelli. Clavton. fc Mathisl l|l989h . We note that this 



be compared with lSavage fc Sembachl l)1996ft who investi- 
gated G/D using HST abundance measurements toward £ 
Oph. When calibrating this measurement with solar abun- 
dances, they found G/D ~ 72. This number corresponds 
to diffuse medium at low extinction (Ay ~ 1.2 mag). Our 
derived values of G/D for p Oph is consistent with this 
value within ~ 10 %. We have howev er used other dust 
grain models, like iDesert et all l(l990l) . which use a dif- 
ferent grain composition and slightly different extinction 
cross-sections, and found G/D ~ 13 0. While this value 
seems to be ruled out by the work of ISavage fc Sembachl 
l)l996h . it shows that the conclusion is rather sensitive to 
the grain composition. Also, the inspection of Fig. IA.1I 
shows that a wide range of values of grain size parameters 
(a m i n , a max , index q) is compatible with a large value of 
Ry. 



5.4. Relationship between metallicity and Ry 

Both explanations ("local" metallicity or larger Ry) for 
a lower Ah,x/Aj ratio in p Oph arc plausible and the 
values of metallicity and Ry quoted above correspond to 
extreme cases. Either grain coagulation has occurred and 
Ry ~ 6 f i!5.2|) . or the standard galactic Ry = 3.1 applies 
but the galactic N-^/Ay relation was derived from X-ray 
measurements using the "old" solar abundances. The "lo- 
cal" metallicity is indeed lower than the "old" solar values 
as shown in Figs^l Therefore, we have looked for an inde- 
pendent indication in p Oph from the CO /H2 ratio derived 
from millimeter-wave measurements. But as explained in 
Appendix B, CO may be depleted onto grains, and above 



Vuong et aL: Determination of the gas-to-dust ratio in dense clouds using X-ray absorption measurements 



17 







[SI 



1.2 



.0 



0.8 



0.6 



1 ' 'I ' 


y< Z"local" 






/_. ; . . . 




# ; / 




i / 
i / 
i ' 
1 / 
' / 
i / 
* / 
/ / 
,» /, 




p Oph 



4 
Rv 



Fig. 11. Metallicity Z as a function of Ry for p Oph 
using the "old" solar and "local" abundances. The 
shaded area represents the uncertainties. The dashed 
line is the extension of the empirical relation of 
ICardelli. Clavton. fc Mathisl ilQSflh valid for 2.6 < Ry < 
5.6. The dotted and dashed-dot ted lines indicate the 
meta llicity for Ry = 3.1 and 4 l|Vrba. Covne. fc Tapial 
Il993h using the "local" and "old" solar abundances. 

all the uncertainties on this ratio are much larger than the 
~ 20—30% accuracy needed to test metallicity values. 

The trade-off between these two interpretations of our 
results - Ry vs. Z - can be expressed in a quantitative 
manner. We now combine them to show they compen- 
sate each other: this will be helpful when more obser- 
vational constrai nts on the metallicity and/q r Ry come 
about. Using the ICardelli. Clavton. fe Mathisl {l989) em- 
pirical relation between Ry and Aj/Ay (found for Ry = 
2.6 to 5.6), we obtain a relation between the cloud metal- 
licity relative to solar, Z/Z@, and Ry : 



Z 

2© 



(JV H Mj)pOph 
(JVkMv)gd 



0.4008 



0.3679 



Ry 



(1) 



taking Nn/Ay = 1.9 x 10 21 cm 2 mag -1 throughout the p 
Oph cloud, and using the fact that the hydrogen column 
density derived from the X-ray spectra is almost inversely 
proportional to the metallicity Z for 0.4 < Z/Z & < 1.2. 
Beyond this range, the higher order terms of the -/Vr-.x — Z 
relation can no longer be neglected. 

This Z(i?v) relation is plotted in Fig. ^2 using the 
"old" solar and "local" abundances. The solid portion 
of the curve represents the most likely region where the 
p Oph dense gas lies, constrained by Z ~ Z\ oca \ and 
Ry = 3.1. The shaded area shows the uncertainties on 
Ry (resulting from the uncertainties on -/Vh,x/^4j, ~ 10%) 
and on Z/Z Q (from Eq.(l)). As discussed qualitatively 
above, the p Oph metallicity is close to the "new" solar 
abundances for Ry close to 3.1 when assuming that the 
galactic gas-to-dust ratio remains unchanged. This is con- 
sistent with the local ISM abundances toward B stars and 
photospheric F & G star abundances. 



6. Conclusions and implications 

We have used the X-ray spectra of 41 young stars without 
circumstellar matter (Class III IR sources, or equivalently 
"weak-line" T Tauri stars) obtained with Chandra and 
XMM-Newton, to derive the hydrogen column densities 
Ar-.x in six nearby molecular clouds, and we have deter- 
mined or re-determined the corresponding dust extinction 
for these stars in the J— band, Aj. Only in the p Oph cloud 
do we have a large enough sample (19 stars) to establish 
a linear correlation between Ah,x and Aj: 

Nn,x/Aj(Z eold ) = 5.6(±0.4) x 10 21 (cm- 2 mag- 1 ), 
JV H ,xMj(Zlocal) = 7.2(±0.5) x 10 21 (cm- 2 mag- 1 ), 

using the "old" solar and "local" abundances for the X- 
ray extinction, respectively. These relations are valid up 
to Aj < 14 (Ay < 50), i.e., probe for the first time the 
gas-to-dust ratio in the dense interstellar medium up to 
very large extinctions. Also, we find that there is no sig- 
nificant deviation of the Nn t x/ Aj ratio along the various 
lines of sight toward the X-ray sources, i.e., no significant 
difference as a function of their respective column densi- 
ties. 

Comparing with the similar relation for the Galaxy 
(Table 5), (A H /^j) g ai= 6.4-7.8 x 10 21 cur 2 mag" 1 , ob- 
tained by using the usual value of the total-to-selective 
extinction ratio Ry = 3.1 (which gives Aj/Ay = 0.28), 
we find that the p Oph AH.xA4j(^©oid) relation lies sig- 
nificantly (> 2a) below the galactic relation. 

We show that this result is consistent with the recent 
downwards revision of the solar abundances. We find that 
the lower value of (jVH/^4j)pOph previously determined us- 
ing the "old" solar abundances changes by 20% when us- 
ing the new set of abundances. But the galactic relation 
derived from the X-ray absorption (see Table EJ) was ob- 
tained using "old" solar abundances. Therefore, the differ- 
ence between the p Oph and Galactic relation persists and 
can be accounted for entirely by a difference in metallicity. 

In order to compare with the galactic relation derived 
from previous X-ray absorption measurements, we keep 
the "old" abundances. We investigated the possibility that 
this was due to a change in Ry, resulting from a change 
in dust grain properties (size distribution, composition, 
etc.). We find Ry = 6.0 ± 2.5. 

Quantitatively, the combination of the dark cloud 
metallicity Z/Zq and of Ry (Fig. 1111) shows that the 
metallicity in p Oph is close to the "local" abundances 
(the Sun, photospheric F & G stars and diffuse ISM to- 
ward B stars) for Ry close to 3.1 if the galactic gas-to-dust 
ratio remains unchanged. 

By contrast, a similar study on five other star- forming 
regions observed in X-rays turned out to have a limited 
value, with only a few X-ray sources with sufficiently high 
count statistics, and insufficiently well-defined IR prop- 
erties. These studies do not go deeper than Aj ~ 1.5 
(Ay ~ 5), preventing us to derive a reliable -/Vh,x/Aj ra- 
tio. 
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The use of the updated solar metallicity and the rela- 
tive difference between the local (d < lkpc) and the galac- 
tic ISM at large distances (up to 4 kpc) gas-to-dust ex- 
tinction relation has some interesting implications: 

1. A first implication is on the chemical evolution of 
the Galaxy. By cross-correlating X-ray absorption and 
abundance measurements, we have shown that the 
metallicities of the Sun, of nearby F & G stars, of 
the local diffuse ISM and dense ISM (in the case 
of the p Oph cloud) are all identical within ~ 10%. 
This seems a priori paradoxical, since the associated 
timescales for the nucleosynthetic evolution of the 
Galaxy are very different, being respectively 5 Gyr, 
~ 1 Gyr, a few Myr, and ~ 1 Myr: we would ex- 
pect the youngest objects (namely the dense ISM) to 
be chemically the most evolved, i.e., to have larger 
metallicities. This "lazy nucleosynthesis" is in con- 
trast with the strong chemical enrichment observed 
over kpc scales towards the center of the Galaxy. 
Mis hurov. Lenine. fc Acharoval l(2002h have noted the 
almost fiat behavior of the metallicity as a function 
of galactic radius in the solar vicinity (which is now 
confirmed by our measurement of the metallicity of 
the p Oph cloud). According to the chemical diffusion 
model proposed by these authors, this can be explained 
by a slow mixing and low star formation rate at the 
galactic solar radius, resulting from to the fact that the 
solar radius lies close to the galactic corotation radius. 

2. A second (related) implication is on the metallicity 
of the Galaxy. We have shown using X-rays that the 
metallicities of the p Oph cloud (using the local abun- 
dances) and of the Galaxy (also using X-rays, but with 
the 'old" solar abundances) differ by > 20%. This 
can be interpreted as the Galaxy being "overmetal- 
lic" with respect to the p Oph cloud and the local 
medium (stars + ISM). As recalled above, the Galaxy 
does show an enhanced metallicity towards the galac- 
tic center. This means that (i) the X-ray absorption 
measurements are indeed quite sensitive to the ISM 
metallicity provided the absorption column density is 
high enough, and (ii) that by a numerical coincidence, 
the ~ 20% downwards revision of the local abundances 
matches the metallicity difference between the local 
medium and the galaxy averaged over the long lines 
of sight of the galactic X-ray absorption measurements 
(see Tables n an d EJ- A detailed modeling of the metal- 
licity along the directions towards distant galactic X- 
ray sources (compact binaries and SNR, u p to ~ 4 
kpc), using iMishurov. Leoine. fc Acharoval {2002) for 
example, would be required to demonstrate this quan- 
titatively, but this is beyond the scope of the present 
paper. 

3. A third implication is on grains in dense clouds. The 
main conclusion is that no strong deviations from the 
usual grain properties are found, like size or compo- 
sition, even in dense regions, as constrained by the 
behavior of Fly using grain models. However, the un- 



certainties are large, and some qualitatively important 
changes (like coagulation in the densest regions, CO 
depletion on grains, etc.) could go unnoticed. But if 
one could independently and precisely measure the 
cloud metallicity (for instance by measuring the pho- 
tospheric metallicity of the young stars they gave birth 
to), one could estimate Ry in dense regions from the 
curve in Fig. ^2 This may then better constrain the 
grain properties. As far as we are aware, there has 
been up to now only one attempt a t measur i ng th e 
metallicity of young low-mass stars bv lPadgettl £996). 
She determined the iron abundances of about 30 PMS 
stars in nearby clouds, and found their [Fe/H] to be 
solar within ±0.15 dex: this uncertainty however is 
larger than the precision level of ~ 10% which we need 
on the metallicity, so more accurate measurements are 
needed. 

4. A fourth implication is more practical, and has to do 
with the determination of X-ray luminosities. In recent 
years, as illustrated in this paper, X-ray spectra from 
the new generation of X-ray satellites have allowed to 
directly measure hydrogen column densities, which in 
turn allow to correct the X-ray luminosities for ex- 
tinction. In all cases of the existing literature, XSPEC 
fits were made by default using "old" solar ISM abun- 
dances. Our work shows that XSPEC should use abun- 
dances revised downwards by ~ 20% for distances up 
to < lkpc, with the resulting (upwards) corrections on 
Lx for sources in nearby clouds. 



The present work should be considered as a pilot study, 
illustrating various successes and problems which can be 
attacked using X-ray absorption in molecular clouds in re- 
lation with IR extinction data. More work, on more clouds, 
including new X-ray and IR observations, for which exist- 
ing archival data are clearly not optimized or insufficient, 
is needed to extract all the potentialities of the method. 
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Appendix A: Modeling the interstellar grain 
population 

A.l. Modeling 

We model N H /A,s and Ry using the lDraine fe Led l)l984f) 
extin ction cross section and the MRN gr ain size distribu- 
tion l|Mathis. Rumpl. fe NordsieckL fl9T7^ . 

iVn/Aj depends on the gas-to-dust mass ratio 
Afgas/Mdust (hereinafter G/D) and the extinction cross 
section Qj iCX t(a) averaged over a size distribution of grains 
n(a), where a is the grain (assumed spherical) radius. 
Absorption and scattering cros s sections are calcu lated 
for grains of radius up to 1 fim l|Praine fc Lee l. ll984h . We 
consider a mixt ure of silicate (53%) and graphite (47%) 
as described bv lDraine fc Leel l)l984h . The grain size dis- 
tribution is described by a power law: n(a)da oc a~ q da, 
<2min < a < flmax, and n(a)da is the density of grains 
along the line of sight with radii in the range a to 
a + da and q =3.5, amin=0-005 urn and a max =0.25 \xm 
(|Mathis. Rumpl. fc NordsieckL Il977j) . Specifically, 



MRN 



M gas 1 
M dust m H 



M, 



dust 



jOmax Q J ext ( a ) 7ra 2 n ( £ j)^ a 



with M, 



dust 



p(a)—7ra 3 n(a)da. 
3 



(A.l) 



(A.2) 



The dust mass density p(a) is assumed to be 3.3 g cm -3 
for silicate and 2.26 g cm -3 for graphite. 

Ry only depends on the extinction cross section at V 
(0.55 pim) and B (0.44 fim) bands respectively, averaged 
over the grain size distribution: 



Ry = 



< a 2 <2v : cxt(a) > 



?B,ext(a) > - < a 2 (3v,oxt(a) > 



(A.3) 



We note that (i) Ry does not depend on the ratio 
G/D and (ii) Ry is very sensitive to the grain composition 
(graphite or silicate) and to the dust size distribution. The 
simplest way to obtain Ry is to measure it in B and V 
bands but this is very difficult for sources with Ay > 4 — 5 
mag. One can then use sub-mm observations to constrain 
the size distribution and/or optical properties. 

The mean grain size is defined as follows: 



< a >- 



J° m " an(a)da 
IaZ7 n ( a ) da 



(A.4) 



The solid curves in Fig. IA.lL b, c show Ry (a), the 
gas-to-dust extinction ratio N-r/Aj (b) and the mean grain 
size (c) as a function of a m j n as derived from Eqs (A.3), 
(A.l) and (A.4), respectively. The two dashed lines in Fig. 
lA.lb represent the values of iVn,x/^4j for Ry — 6 (or 
Z/Zq = 1) and Ry = 4 (or Z/Z Q = 0.9) (see Fig.LTTJ. 
For Ry between 4 and 6, we find that G/D is comprised 
between 95 and 80 and a m ; n between 0.025 and 0.07 /mi. 
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Fig. A.l. Ry (a), gas-to-dust extinction ratio Nn/Aj (b), 
and the mean grain size (c) as a function of a m i n . The 
dashed line and the shaded area in (b) are the mean value 
±lcr of the observed -/Vh,x/^4j in P Oph. The dotted lines 
represent the interval of a m ; n for Ry between 4 and 6. 



Appendix B: The N(CO)/N(H 2 ) ratio and CO 
depletion onto grains 

We discuss here whether we can test the conclusion from 
our combined X-ray/IR analysis that the ISM abun- 
dance of heavy atoms (or metallicity) is lower than so- 
lar in p Oph, using the CO/H 2 ratio derived from X- 
ray and millimeter observations. We take as a test case 
WL 5, the most deeply extincted source in our sample 
(Ay ~ 45; see §4.2.1), using the total hydrogen (mostly 
H 2 ) c olumn density Nu y de rived from the X-ray spec- 
trum. IWilking fc Ladal Il983) have done C 18 observa- 
tions of the dense cores of the p Oph cloud. They de- 
termined the column density of C 18 adopting a local 
thermodynamic equilibrium (LTE) approach. Along the 
line of sight toward WL5, the C 18 Q column density is 
1.53xl0 16 mag.cm" 2 llWilking fc Ladal Il983l) . Since H 2 
is the dominant component of H along the line of sight, 
the H 2 column density is N(H 2 )=A^h,x/2. To convert 
N(C ls O) to N( 13 CO), we use the isotop ic abundance ra- 
tios J 12 C/ 13 C1=77 and [ 16 O/ 18 O]=560 llWilson fc Roodl 
Il994h . For WL5, we find N(H 2 )/N( 13 CO)~3.3 xlO 5 us- 
ing solar abundances, and N(H 2 )/N( 13 CO)~4.5 xlO 5 us- 
ing B-star abundances. These values are consistent with 
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the usual relatio nship N(H 2 )/N( 13 CO)=5.0±2.5 x 10 5 by 
iDickmanl l)l978|) . Note that the determination of N(CO) 
does not depend on metallicity. In conclusion, given the 
errors we are not sensitive to changes in the CO/H 2 ra- 
tio, whether we take solar or B-star abundances. Since 
we are looking for evidence of underabundances, the fact 
that CO can be depleted o n grains in dense clouds (see 
lAlves. Lada. fe Ladal <|l999j) for a study of the correlation 
between C 18 line emission and dust extinction toward 
the L 977 molecular cloud) is an additional difficulty. We 
can only exclude an underabundance of CO (intrinsic or 
by depletion) by a factor > 2, i.e., much larger than the 
effect we are looking for. 
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